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[Method for Magnetic Resonance 
Fluid Characterization] 

Background of Invention 

[0001] The Invention relates generally to the field of well logging. More particularly, the 
invention relates to improved techniques for well logging using nuclear magnetic 
resonance and methods for analyzing the nuclear magnetic measurements. 

[0002] Oil well logging tools include nuclear magnetic resonance (NMR) instruments. NMR 
instruments can provide a wealth of information for formation evaluation that is not 
obtainable from other well logging measurements. Information provided by NMR 
measurements include the fractional volume of pore space, the fractional volume of 
mobile fluid filling the pore space, and the porosity of earth formations. General 
background of NMR well logging is described in U.S. Patent No. 6.1 40,81 7A1, 
assigned to the assignee hereof. 

[0003] 

The signals measured by nuclear magnetic resonance (NMR) logging tools typically 
arise from the selected nuclei present in the probed volume. Because hydrogen nuclei 
are the most abundant and easily detectable, most NMR logging tools are tuned to 
detect hydrogen resonance signals (form either water or hydrocarbons). These 
hydrogen nuclei have different dynamic properties (e.g., diffusion rate and rotation 
rate) that are dependent on their environments. The different dynamic properties of 
these nuclei manifest themselves in different nuclear spin relaxation times (i.e., spin- 
lattice relaxation time ( T ^ ) and spin-spin relaxation time ( )) and diffusion 
constants. For example, hydrogen nuclei In viscous oils have relatively short relaxation 
times and low diffusivity, whereas hydrogen nuclei in light oils possess relatively long 
relaxation times and high diffusivity. Furthermore, the hydrogen nuclei In free water 
typically have longer relaxation times than those in bound water. Consequently, these 
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differing NMR relaxation times can provide information on properties of the earth 
formations. 

[0004] Most NMR logging tools measure the spin-spin relaxation times ( 7"^ ) to derive 
the properties of the earth formations. relaxation Is often measured from a train 
of spin-echoes that are generated with a series of pulses such as the Carr-Purcell- 
Meiboom-Gill (CPMC) pulse sequence or some variants of this. The CPMG pulse 
sequence is well known in the art. See Meiboom, S., Gill, D., 1958, "Modified Spin Echo 
Method for Measuring Nuclear Relaxation Times, "Review of Scientific Instruments, 29, 
688-91. 

[0005] Although measurements provide useful information for formation 
characterization, Information alone may not be sufficient for distinct 
characterization of earth formations, especially when different components In the 
formations have similar or overlapping values. In order to differentiate different 
fluids (e.g., hydrocarbons versus connate water) with similar or overlapping 
distributions, several differential methods have been proposed, such as the 
differential spectrum method (DSM) and time domain analysis (TDA). These methods 
takes advantages of different longitudinal relaxation times of different fluids. Two 
sets of measurements are made with different wait times (times for the spin to be 
polarized by the static magnetic field). One wait time is selected that one type of fluid 
(e.g., brine that has a longer longitudinal relaxation time ) would not fully relax. 
As a result, the signal magnitudes from the long fluid would be substantially 
reduced. Subtraction of these two sets of measurements would then produce a 
difference measurement which is comprised mostly of signals from the fluids with 
long times. However, success of these approaches relies on the selection of proper 
wait times, which requires prior knowledge of the NMR properties of the fluids in the 
formation. 

[0006] 

Recently, a magnetic resonance fluid characterization (MRF) method has been 
shown to provide more useful information. For a detailed discussion of the MRF 
method, see U.S. Patent No. 6,229,308 Bl issued to Freedman and assigned to the 
assignee of the present invention. This patent is hereby incorporated by reference. 
When T distributions overlap, the MRF method distinguishes oil and water in porous 
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media based on different molecular diffusion. For the same , oil and water have 
different diffusion constants. Therefore, contributions of different fluids to the 
measured 7"^ distributions can be separated by combining measurements with 
different sensitivity to diffusion. However, in order to use MRF, diffusion constants 
must be measurable. If becomes comparable with or shorter than the decay due to 
diffusion, then the diffusion constant can no longer be determined, and the MRF 
method becomes impracticable. In other words, the MRF method is applicable only 
when values are longer than a certain critical value. 

[0007] While the MRF analysis has proved to be a powerful approach, it is desirable to 
have methods that can be used to analyze fluids with not only long T ^ or T ^ , but 
also short ^ ^ ^2 ' ^"^^ viscous fluids. Furthermore, it is desirable to have 
methods that can be used in logging tools with a low magnetic field gradient or a 
saddle point in the magnetic field. 

Summary of Invention 

[0008] One aspect of the invention relates to methods for determining properties of a 
mixture of fluids using NMR data that include longitudinal and transverse relaxation 
information. According to embodiments of the invention, a method for determining 
properties of a mixture of fluids includes: (a) acquiring a plurality of nuclear magnetic 
resonance measurements from the mixture of fluids, each of the plurality of nuclear 
magnetic resonance measurements having a different value in an acquisition 
parameter for which at least one relaxation selected from the group consisting of 
longitudinal relaxation and transverse relaxation affects magnitudes of the nuclear 
magnetic resonance measurements; (b) generating a model of the mixture of fluids; 
(c) calculating a synthesized nuclear magnetic data set based on the model; (d) 
comparing the synthesized nuclear magnetic data set with the nuclear magnetic 
resonance measurements; and (e) adjusting the model and repeating (c) and (d), if 
difference between the synthesized nuclear magnetic data set and the nuclear 
magnetic measurements is greater than a minimum. 

[0009] 

Another aspect of the invention relates to methods for logging an earth formation 
surrounding a wellbore. According to embodiments of the invention, a method for 
logging an earth formation surrounding a wellbore includes: (a) lowering a nuclear 
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magnetic resonance instrument into the wellbore; (b) inducing a static magnetic field 
in a region of investigation; (c) generating a series of radio frequency magnetic field 
pulses in the region of investigation, and receiving signals comprising a train of 
nuclear magnetic resonance spin echoes in response to the series of radio frequency 
magnetic field pulses, wherein the generating and the receiving are repeated a 
plurality of times each with a different value in an acquisition parameter for which at 
least one of the longitudinal relaxation and transverse relaxation affects magnitudes 
of the signals; (d) generating a formation model that includes at least one component 
for a connate water phase and at least one component for an oil phase; (e) calculating 
a synthesized nuclear magnetic data set based on the formation model; (f) comparing 
the synthesized nuclear magnetic data set with the nuclear magnetic resonance 
measurements; and (h) adjusting the formation model and repeating (c) and (d). If 
difference between the synthesized nuclear magnetic data set and the nuclear 
magnetic measurements Is greater than a minimum. 

[0010] Other aspects and advantages of the invention will be apparent from the following 
description and the appended claims. 

Brief Description of Drawings 

[001 1] FIG. 1 is a diagram of a typical setup for well logging. 
[001 2] FIG. 2 shows a diagram of a prior art NMR logging tool. 

[001 3] FIG. 3 is a block diagram of circuitry for producing the RF pulses and receiving the 
spin echoes. 

[0014] FIGs. 4A and 4B show the exponential recovery of signals after an inversion pulse 
or a saturation pulse. 

[001 5] FIGs. 5A and SB show the pulse sequences for IR-CPMG and SR-CPMG. 
respectively. 

[0016] FIG. 6 shows a flow chart of a method according to one embodiment of the 
invention. 

[0017] 

FIGs. 7A and 78 show the results of computed porosity and T , respectively, 
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versus connate water T ^ /r^ values. 

[001 8] FIG. 8 shows a graph of quality of model fit to experimental data as a function of 

T /T ratios. 
1 2 

Detailed Description 

[0019] FIG. 1 shows a schematic of a nuclear magnetic resonance (NMR) logging tool 30 
for investigating earth formations 31 traversed by a borehole 32 . The NMR logging 
tool J^7is suspended in the borehole 32 ox\ an armored cable 33 , the length of which 
substantially determines the relative depth of the logging tool 30 . The cable length is 
controlled by suitable means at the surface such as a drum and winch mechanism 8 . 
Surface equipment Zmay include a processor subsystem to communicate with 
downhole equipment Including NMR logging tool 30 . 

[0020] The NMR logging tool Jf7can be any suitable nuclear magnetic resonance logging 
device. It may be one adapted for in wireline logging applications as shown in FIG. 1 , 
or one that can be used in logging while drilling (LWD) applications. The NMR logging 
tool J^?typically includes a means for producing a static magnetic field in the 
formations, at least one radio frequency (RF) antenna, and means for producing pulses 
of RF power to induce RF magnetic fields in the formations and for receiving the spin 
echoes from the formations. The means for producing a static magnetic field may 
comprise a permanent magnet or magnet array, and the RF antenna may comprise one 
or more RF antennas, which may be solenoid antennas, loop antennas, or saddle 
antennas. 

[0021] 

A schematic representation of some of the components of an NMR logging tool 30 
is illustrated in FIG. 2, which shows a first centralized magnet or magnet array 36 and 
an RF antenna 37 , which may be a suitably oriented coil or coils. FIG. 2 also illustrates 
a general representation of closely-spaced cylindrical thin shells, 38-1 , 38-2 . . . 38- 
A/, that can be frequency selected in a multi-frequency logging operation. These thin 
shells are resonance regions where magnet J6^ has a radial gradient in the field 
amplitude. One such device is disclosed in U.S. Patent No. 4,710,713 issued to 
Strikman. In FIG. 2, another magnet or magnet array 39 \s shown. Magnet array 39 
may be used to pre-polarize the earth formation ahead of the investigation region as 
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the logging device 30 \s raised in the borehole in the direction of arrow Z. Examples of 
NMR logging tools having pre-polarization magnet arrays are disclosed in U.S. Patent 
Nos. 5,055,788 and 3,597,681 . 

[0022] FIG. 3 illustrates a schematic of a circuitry of an NMR tool for producing the RF 
pulses and for receiving and storing the spin echoes. One skilled in the art would 
appreciate that any other suitable circuitry could be used without departing from the 
scope of the invention. In FIG. 3, a downhole processor 270has associated memory, 
timing, interfaces, and peripherals (not separately shown), as known in the art. The 
processor subsystem 210 is coupled with telemetry circuitry 20S , for communication 
with a processor (not shown) at the earth's surface. The pulse forming circuitry 
includes a variable frequency oscillator 220vifh\ch, under control of processor 210 , 
produces radio frequency (RF) signals at the desired frequencies. The output of 
oscillator 220 is coupled to a phase shifter 222 , which permits control of pulse 
phases, and then to a modulator 230 , both of which are under the control of 
processor subsystem 270 . The phase shifter 222 and modulator 230 ca.n be 
controlled, in a manner known in the art, to produce the desired pulse phases of RF 
field. The output of modulator 230 is coupled, via a power amplifier 23S , to the RF 
antenna 240 . A Q-switch 2S0 can be provided to damp the RF antenna system to 
reduce antenna ringing. The antenna 240 \s also coupled with a receiver section via 
duplexer 26S , the output of which is coupled to receiver amplifier 270 . The duplexer 
265 protects the receiver amplifier 270 from the high power pulses which pass to the 
RF antenna 240 during the transmitting and damping modes. During the receiving 
mode, the duplexer 25S acts as a low impedance connection from antenna 240 to the 
receiver amplifier 270 . The output of receiver amplifier 270 is coupled to a dual 
phase-sensitive detector 27S , which also receives, as a reference, a signal derived 
from the oscillator signal. The detected output is coupled to analog-to-digital 
converter 280 , the output of which is a digital version of the received nuclear 
magnetic resonance signal. Although the logging device or tool 30 \s shown as a 
single body In FIG. 1, it may alternatively comprise separate components, and the tool 
may be combinable with other logging tools. Also, while a wireline is illustrated, 
alternative forms of physical support and communicating link can be used, for 
example in a measurement while drilling system. 
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[0023] Several NMR parameters may be measured that can be used to derive formation 
properties. Most NMR logging operations measure the spin-lattice (longitudinal) 
relaxation times ( 7"^ ) and/or spin-spin (transverse) relaxation times ( 7"^ ) of 
hydrogen nuclei. In addition, some NMR logging tools may provide a ratio of T ^ / T ^ 
directly, and other NMR tools may provide diffusion constants ( D). 

[0024] Various pulse sequences are known in the art for measuring the NMR relaxation 
times. For example, T ^ (longitudinal) relaxation may be measured using an 
inversion-recovery (IR), saturation recovery (SR), or steady-state approach, or any of 
their derivatives. T ^ relaxation is often measured from a train of spin-echoes that are 
generated with a series of pulses such as the Carr-Purcell-Meiboom-Gill (CPMG) pulse 
sequence or some variant of this. The CPMG pulse sequence is well known in the art. 
See Meiboom, S., Gill, D., 1958, "Modified Spin Echo Method for Measuring Nuclear 
Relaxation Times," Review of Scientific Instruments, 29, 688-91. CPMG pulse 
sequence generates a train of spin echoes, whose amplitudes exponentially decay as a 
function of time. The exponential decay life time is referred to as a transverse (spin- 
spin) relaxation time, 7"^ . Thus, T ^ measurements are accomplished by analyzing 
the amplitudes of spin echoes thus obtained. 

[0025] 

NMR measurements of diffusion constants are accomplished in the presence of 

magnetic field gradients. Magnetic field gradients produce different strengths of 

magnetic field at different locations. The different magnetic field strengths manifest 

themselves as different Larmor frequencies in the detected signals (because o) = y B 

o 

). If a delay time is inserted between the pulse and the start of signal collection, 

o 

nuclear spins can diffuse from location A to location B, i.e., from one Larmor 
frequency to another, with a concomitant loss of phase coherence. This results in a 
decrease in the magnitudes of the detected signals. The diffusion induced decay of 
the NMR signal depends on the square of the echo spacing and on the magnetic field 
gradient, because the decay rate is governed by: 

where D is the diffusion constant, G is the gradient of the static magnetic field, y is 

the gyromagnetic ratio, and t is the echo spacing. See, e.g C.P. Slichter, Principles 

e 

of Magnetic Resonance, 3 Edition, Springer Verlag, Berlin, Heidelberg, New York, 
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1 990, chapter 8 and appendix G. 

[0026] Once NMR data, which include information on , relaxation times, 7"^ / 

ratio, or diffusion rates D, are collected, they are analyzed with an inversion method 
to derive a model of the earth formations. For the derived model, properties of the 
formation then can be calculated. Many methods are known in^the art for inverting 
NMR data, including the Window Processing (WP) method disclosed in U.S. Patent No. 
5,291 ,1 37 issued to Freedman. However, when different fluids have similar NMR 
parameters, for example, same or overlapping distributions, most inversion 
methods cannot separately characterize these fluids. 

[0027] Freedman et al. disclosed a magnetic resonance fluid (MRF) characterization 

method that is capable of distinguishing different fluids even if they have overlapping 
NMR parameters (e.g., overlapping distributions). See Freedman et aL, " A New 
Method of Fluid Characterization In Reservoir Rocks: Experimental Confirmation and 
Simulation Results , "Paper SPE 6321 4 presented at the 2000 Society of Petroleum 
Engineers Annual Technical Conference and Exhibition, and Freedman et aL, " Field 
Applications of a New Nuclear Magnetic Resonance Fluid Characterization Method 
paper SPE 71 71 3 presented at the 2001 Society of Petroleum Engineers Annual 
Technical Conference and Exhibition. In addition, U.S. Patent 6,229,308 Bl issued to 
Freedman discloses the application of this method in formation evaluation. 

[0028] The MRF method is based on two key concepts: (1) a new microscopic CVM 
(constituent viscosity model) that relates NMR relaxation times and molecular 
diffusion coefficients in crude oils, and (2) a new multi-fluid relaxation model. The 
MRF method provides a detailed formation evaluation of the near wellbore region 
investigated by modern NMR logging tools. The information provided by MRF includes 
flushed-zone fluid saturations and volumes, total and bound-fluid porosities, bulk 
volumes of hydrocarbon saturations, oil viscosities, and hydrocarbon-corrected 
permeabilities. 

[0029] CVM relates individual constituent diffusion-free relaxation times and diffusion 
coefficients to a distribution of constituent viscosities. The constituent viscosities are 
molecular variables that are analogous to the "friction coefficients" used in Langevin 
equation models of Brownian motion in viscous media. Before looking at the 
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correlation between proton relaxation times and diffusion constants and the viscosity 
of a mixture, it is helpful to look at the simplest situation first, i.e., a pure liquid. 

[0030] In pure liquids, Bloembergen et al. in " Relaxation Effects in Nucfear Magnetic 

Resonance Absorption , "Physical Review (1948), 73, no. 7, pp. 679-712, showed that 
proton relaxation times (e.g., T ^ ) depend on the temperature and the liquid viscosity 
according to: 

where a is a constant, T ^ is the spin-spin (transverse) relaxation time, T\s the 

temperature of the liquid in degrees Kelvin, and r| is the viscosity of the liquid. This 

o 

relationship is also valid for T , provided that the "fast motion" condition, oj < 

1 O T 

1 , is satisfied, where t is a rotational (i.e., molecular tumbling) correlation time and o) 
is the Larmor frequency. In pure liquids, the relationship between diffusion 

o 

constants and viscosities is described by the Einstein-Stokes model, which predicts 
that diffusivity is related to temperature and viscosity according to the equation: D = 
/cT/6 TT n i^ ^ where Ic is the Boltzmann constant, R is the radius of the spherical 
particle, n is the viscosity of the fluid, and Fis the temperature in degrees Kelvin. This 
relationship may be simplified as D = bT/ n , where b = k/6 tt R. 



[0031] 



The dynamic properties of each individual constituent in a mixture are similar to 

those of a pure liquid. Thus, a similar relationship exists between the relaxation time 

( T ) and the diffusion constant ( D ) of the Xr-f/? constituent in a mixture of 
2o,k o,k 

oils according to CVM. The subscripts " denotes the parameter is for the Ic-th 
constituent in the oils. The key CVM equations for a hydrocarbon mixture with n 
components are: 

bT 

X?.^=— (2) 
71. =%^*l7/'*-*^„^ ^{vX (3) 

where T . and D are the relaxation time and the diffusion coefficient, 
2o,k o,k 

respectively, for the ic-t/i molecular constituent in a hydrocarbon mixture (e.g., a 
crude oil), T\s the sample temperature in degrees Kelvin, f(GOR) , as will be explained 
below, is an empirically determined "universal function" of gas/oil ratio (COR) and is > 
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1 , a and b are constants, /? is the constituent viscosity for the k-t/i molecular 

k 

constituent, f Is the proton fraction of each of the mixture constituents, and /y Is 
k o 

the macroscopic viscosity ( rj ) of the mixture. Equation (3) indicates that the 

o 

macroscopic viscosity of the mixture equals the logarithmic mean (or geometric mean) 

of the constituent viscosity distribution (( /; , 

k LM 

[0032] The factor f(GOR) is included in equation (1) because it has been shown that GOR 
is an important parameter in determining the relaxation time dependence on viscosity 
and temperature. See Lo et a!., Relaxation Time And Diffusion Measurements of 
Metftane And N-Decane Mixtures , The Log Analyst, pp. 43-46, November-December, 
1 998; see also U.S. Patent No. 6,229,308 Bl . Note that the f(GOR) function is always 
equal or greater than 1. For dead oils (which are oils having no solution gas) and 
mixtures that do not contain any solution gas, f(COR) = f(0) = 1 . (Freedman, R., et al., 
" A New Method of Fluid Characterization in Reservoir Rocks: Experimental 
Confirmation and Simulation Results , "paper SPE 6321 4 presented at the 2000 SPE 
Annual Technical Conference and Exhibition, Dallas, 1-4 October; hereafter referred 
to as "Freedman et al."). The constant a in equation (1) has been empirically 

determined to be about 0.004 s • cp • K ^ and this value has been shown to be 
applicable to substantially every constituent in crude oils. See Freedman et al. The 

-8 

constant b \x\ equation (2) has been empirically determined to be about 5.05 x 10 

cm ^ • s ^ • cp • K V This value for constant b is also universally applicable to 
constituents in both dead and live hydrocarbon mixtures and crude oil. See Freedman 

et al. 

[0033] 

Combining equations (1) and (2) provides the expression: 

Equation (4) indicates that there is a linear relationship between the diffusion constant 

( D ) and relaxation time ( T ) for each of the constituents in the mixture. This 
o,k 2o,k 

linear relationship is similar to the relationship described above for the pure liquids. 

Furthermore, CVM should converge with a pure liquid model when there is only one 

constituent In the mixture (i.e., pure liquid, n = Ar = l). Indeed, when n =1, 

o o 

equations (1) and (2) reduce to a form identical to that for pure liquids as described by 
Bloembergen et al. 
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and Stokes Einstein model 

[0034] The validity of equations (1) - (4) has been shown (as disclosed in Freedman et al.) 
using experimental measurements of T ^ and Z7 distributions on live and dead 
hydrocarbon mixtures and crude oils. Equations (1) and (2) provide some Interesting 
insights into the nature of the molecular interactions in hydrocarbon mixtures such as 
crude oils. First, each molecular constituent retains its separate properties in the 
mixture. Second, the relationships between NMR relaxation times, molecular diffusion 
coefficients, and constituent viscosities for each constituent in a mixture are the same 
as that for a pure liquid, except that the constituent viscosity rather than the pure 
fluid viscosity (macro viscosity) determines the relaxation time and the diffusion 
coefficient of each constituent. Third, all unknown factors that might affect molecular 
interactions, e.g., molecular composition and molecular sizes, are contained in the 
constituent viscosity terms in this model, and there is no need to determine these 
factors. In addition, equation (4) predicts that the T ^ and Z7 distributions in 
hydrocarbon mixtures and crude oils are not independent. 

[0035] The importance of the CVM for fluid characterization can be appreciated from 

equation (4), which correlates diffusion-free relaxation times with molecular diffusion 
in crude oils. This link reduces the number of unknown parameters in the NMR 
multifluid relaxation model and results in robust and accurate recovery of oil T ^ and 
Z7 distributions using the MRF method, even when the oil and connate water (or brine) 
T ^ distributions are broad and largely overlap one another. 

[0036] The MRF method uses a general relaxation model for a formation containing brine, 
hydrocarbons, and oil-based mud filtrate (OBMF). Mud having hydrocarbons as the 
continuous phase is commonly used in drilling the borehole and it invades porous 
formations as OBMF. In the absence of magnetic field gradient or diffusion, the 
amplitudes of the spin echoes in NMR measurements decay exponentially as a 
function of time ( T ^ or T ^ ). In the presence of a magnetic field gradient, molecular 
diffusion also contributes to the reduction of the measured NMR signal amplitudes. 
Let 

^; 
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represent the amplitude of the j-th echo acquired during measurement p . Consider 
the following general relaxation model, 



where the first, second, and third terms are connate water, native hydrocarbon, and 
OBMF signals, respectively. This three-phase model explicitly accounts for 
contributions from individual constituents in the connate water and hydrocarbon 
phases, but only assumes an average relaxation time distribution in the OBMF phase. 
The simplified term for OBMF is Justified because experimental measurements in 
OBMF have shown that NMR relaxation time distributions for individual constituents in 
typical OBMF are very narrow and can be described by a single decay exponent. If such 
simplification is not desirable, one may include a summation function, like that for 
connate water or native hydrocarbon, for the OBMF. In addition, one skilled in the art 
would appreciate that this general relaxation model may be modified, for example, to 
remove the OBMF term, if justified (e.g., investigation of core samples removed from 
formations, or LWD-NMR logs in which mud infiltration is negligible), or to include gas 
in the fluids as disclosed in Appendix A of U.S. Patent No. 6,229,308 Bl. 

[0037] The apparent transverse (spin-spin) relaxations ( ) in any of the three phases 

modeled in equation (5) include actual transverse relaxation and signal decay due to 

self diffusion of molecules in the static magnetic field gradient C . For unrestricted 

P 

diffusion in a uniform magnetic field gradient in the connate water, the apparent 
transverse relaxation rates can be written as. 

Here, T in the first term on the right-hand side are a set of relaxation times that 
2,1 

represent the sum of surface and bulk relaxation of the connate water phase. The 
second term is the contribution to the relaxation rate from diffusion, where 

is the proton gyromagnetic ratio and (T) is the temperature dependent self 
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diffusion coefficient of water in units of cm /s. Note tliat equation (6) assumes an 

unrestricted diffusion and a uniform magnetic field gradient C . One skilled in the 

P 

art will appreciate that corrections to D (T) for the effects of restricted diffusion and 

w 

to G for the effects of internal rock gradients can be applied if appropriate, such as 
P 

when the formations have minerals having high magnetic susceptibility. 
[0038] Similarly, the apparent transverse relaxation rates in the native oil 

can be written in the form. 



[0039] 



CD 



where T ( n ) \s the bulk relaxation time associated with amplitude b , in the 
2.0 ' k k 

hydrocarbon relaxation time distribution, and D (/j ) \s a viscosity dependent 

o k 

diffusion constant. 



The measured viscosity ( rj ) reflects a macroscopic transport property of the 

o 

crude oil that determines its flow properties and is the quantity that is used in 

hydrodynamic transport equations like the Navier- Stokes equation. Morriss eta/. 

showed that, for a suite of dead (i.e., not containing dissolved solution gas) crude oils, 

there exists a strong correlation between the logarithmic mean relaxation times of 

their constituents and the measured viscosities. See Morriss et al., " Hydrocarbon 

Saturation And Viscosity Estimation From NI\/fR Logging in The Belridge Diatomite , 

"Paper C presented at the 3 5th Annual Meeting Of The Society Of Professional Well 

Logging Analysis, 1 994. The macroscopic viscosity ( ri ) of live crude oils is 

o 

empirically related to the logarithmic mean 

of the transverse relaxation time distributions by an equation of the form, 

^ T 

where a is an empirically determined constituent constant that has been determined 
by Looyestijn (" Determination of Oii Saturation From Diffusion NMR Logs " Paper SS 
th 

presented at the 37 Annual meeting of the Society of Professional Well Log 
Analysts, 1 996) to be around 250 (i.e., a ^ 250 Ks cp ), for 

in seconds and rj in centipoise and Tthe temperature in degrees Kelvin. Thus, c = 
o 
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for live oil. The empirically derived function accounts for live oils (those 

containing dissolved solution gas) and has been discussed by Freedman et al . in 
paper number 632 1 4 entitled "A New NMR Method of Fluid Characterization in 
Reservoir Rocl<s: Experimental Confirmation and Simulation Results" pr&sQnted at the 
2000 Society of Petroleum Engineers Annual Technical Conference and Exhibition 
meeting. 

[0040] The n ^ terms in equation (7) are microscopic viscosities that reflect the complex 
composition of crude oils. Analogously with equation (8), the constituent viscosities 
are assumed to be related to the components in the relaxation time distribution via 
the same equation, 

(SO 

[0041] The dependence of the relaxation times on viscosity and temperature in equations 
(8) and (9) is consistent with the experimental observations and theoretical 
predictions of Bloembergen, Purcell, and Pound, Relaxation Effects in Nuclear 
Magnetic Resonance Absorption , Physical Review, vol. 73, no. 7, pp. 679-712, 1948. 



[0042] Similar to the Stokes-Einstein equation, the self-diffusion constants for the crude 
oils, 
oils. 



oils, D , are assumed to have the same dependence on F/ n , • Therefore, for crude 
o k 



00) 

2 

where Z?is a constant, D is the measured crude oil diffusion constant in cm /s, and 

o 

T\s the temperature in degrees Kelvin. The empirical constant for crude oils, 

, is given by Freedman et aL In paper number 63214 entitled "A New NMR Method of 
Fluid Characterization In Reservoir Rocks: Experimental Confirmation and Simulation 
^e5^/f5, "presented at the 2000 Society of Petroleum Engineers Annual Technical 
Conference and Exhibition meeting. 

[0043] Analogously to the relationship between the macroscopic diffusion constant ( D ) 

o 
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and the macroscopic viscosity ( q ), the microscopic constituent diffusion constants 

o 

( D ( ri . )) of crude oil mixtures are related to the microscopic constituent viscosities 
o k 

( /; ) according to the following equation: 

K 

(11) 

Equation (1 1) implies that there exists at the molecular level a distribution of diffusion 
constants in the crude oil mixture analogous to the distribution of relaxation times. 
This is the key concept of the CVM (constituent viscosity model). These distributions 
of diffusion constants and relaxation times can be derived from the measured NMR 
data by iteratively fitting a multifluid model as shown in equation (5) to these data 
using the method disclosed in U.S. Patent No. 6,229,308 Bl issued to Freedman. 

[0044] The MRF method inverts suites of NMR data using a multifluid relaxation model as 
shown in equation (5). In this model, different fluids that have overlapping T ^ 
distributions are distinguished based on differing diffusion constants. This is feasible 
because water molecules typically diffuse much faster than oil molecules, and much 
slower than gas molecules. Consequently, NMR diffusion measurements can provide a 
means for detection and differentiation of reservoir fluids, even if other NMR 
parameters cannot. 

[0045] jq the contributions of differing molecular diffusion rates to resolve the fluid 
composition where there are overlapping T ^ distributions, the MRF method uses a 
suite of spin-echo measurements which have been acquired with different echo 
spacings, polarization times, applied magnetic field gradients, and/or numbers of 
echoes. The suite of NMR data may be acquired with a conventional CPMG-pulse 
sequence or a new type of "diffusion-editing"("DE") spin-echo sequence tailored for 
fluid typing. 5eeM.D. H u rlimann et al., " Diffusion-Editing: New NMR Measurement 
of Saturation and Pore Geometry paper presented at the 2002 Annual Meeting of the 
Society of Professional Well Log Analysts, Osio, Japan, June 2-5; see also , U.S. 
Application Serial No. 09/723,803 filed on November 28, 2000 by H u rlimann, 
entitled " Simpler and More Robust Hydrocarbon Typing with Pulsed NMR ." This 
application is assigned to the assignee of the present invention and is hereby 
incorporated by reference. DE sequences are similar to the CPMG sequences except 
that the initial two echoes are acquired with longer echo spacings and the third and 
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subsequent echoes are acquired with shorter echo spacings. Diffusion information is 
encoded during the acquisition of the first two echoes, whereas the third and 
subsequent echoes can provide bulk and surface relaxation time information with 
little, if any, attenuation of the signal by diffusion. 

[0046] The suite of spin-echo measurements are then used in inversion analysis to fit the 
general multifluid relaxation model as shown in equation (5). Any inversion method 
known in the art may be used, including the Window Processing (WP) method 
disclosed in U.S. Patent No. 5,291 ,1 37 issued to Freedman. This patent is assigned to 
the assignee of the present invention and is hereby incorporated by reference. For 
more detailed discussion of the MRF methods see Freedman et al., " A New Method of 
fluid Characterization in Reservoir Rocks: Experimental Confirmation and Simulation 
Results " paper SPE 63214 presented at the 2000 SPE Annual Technical Conference 
and Exhibition, and Freedman et al., " Field Applications of a New Nuclear /[Magnetic 
Resonance Fluid Characterization Method paper SPE 71713 presented at the 2001 
SPE Annual Technical Conference and Exhibition. 

[0047] As stated above, for the diffusion constants to be measurable, the nuclear spin 
relaxing processes (longitudinal or transverse relaxation) should be slower than the 
diffusion process. Otherwise, the NMR signals would have diminished too much before 
any appreciable diffusion has occurred. In other words, if the NMR signal decays too 
quickly such that the decay becomes comparable with or faster than the signal decay 
due to the diffusion process, the diffusion constants can no longer be determined. 
Without the diffusion constants, the fluids with overlapping T ^ distributions cannot 
be separated by the conventional MRF method. 

[0048] 

The lower limits of T ^ or T ^ , above which the MRF method is useful, are even 

TM 

higher for logging-while-drilling (LWD-NMR ) tools because these tools typically 
have low magnetic field gradients. As is evident from the above discussion, the 

magnitudes of the diffusion-induced relaxation depend on the strengths of the 

TM 

magnetic field gradients. With low-gradient NMR tools (e.g., LWD-NMR or MRWD 
TM 

), it requires longer delay times (e.g., the Inter-echo delays) for the diffusion 
information to become detectable in the measurements. The required longer delay 

time makes the diffusion measurement only possible for fluids with relatively long T 
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or T ^ times. For example, a typical LWD-NMR tool operating at a gradient of 4 C/cm 
(Gauss per centimeter) would need about 6 times longer delays, as compared with a 
TM 

typical MRX tool operating at a gradient of 2 5 G/cm, for the molecular diffusion to 
be measurable. Thus, using a low-gradient NMR tool, the diffusion-based MRF 

method is useful only for oils with long ^ -^^^ ^ 2 ^^^^ fluids, but not for 

heavy oils. 

[0049] Embodiments of the present invention are capable of extending the applicability 
of the MRF method to situations where short NMR relaxation times do not permit 
accurate determination of the diffusion constants. Embodiments of the invention use 
other characteristics of the fluids, for example, different ^ ^ I ^ ^ "^^^^^^ ^ ^ '"^ 
equation (5)), to differentiate fluids having similar or overlapping T ^ distributions. T 
^1^2 ^^^^^^ readily measurable using typical NMR logging tools. T ^ / 
contrasts (different T ^ / T ^ ratios between different fluids) have been used in the art 
for other purposes, such as "polarization correction," which is important for 
measurements acquired using pulse sequences that have short wait times (WT). See 
U.S. Patent No. 6,229,308 Bl, issued to Freedman. 

[0050] In free fluids, T ^ equals T ^ except for very viscous fluids (e.g., tar) and/or for 

very high Larmor frequencies (usually not encountered in logging tools). Thus, in free 
fluids, the T ^ / T ^ ratio equals 1. The porous spaces in the formations often are 
water wet. As a result, oils (hydrocarbons) are repelled from the surfaces (walls) of the 
pores, and they behave like free fluids without boundary effects. In contrast, water 
interacts with the surfaces of the rock grains. Consequently, each water molecule 
would spend some time as a surface-bound water molecule and other time as a free- 
water molecule. Surface-bound water molecules tend to have faster transverse (spin- 
spin) relaxation, i.e., short T ^ , while free-water molecules would have longer T ^ 
values that approach the T ^ values. Therefore, on average, water molecules in the 
formations have shorter (average) T ^ values than those of free fluids, e.g., oils 
(hydrocarbons). Consequently, for fluids in the formations, there is a / T ^ 
contrast between oil and water: T ^ / T ^ (oil) ^ 1, and T ^ / (water) > 1, e.g., 1,5 
3.0 in sandstones. 



[0051] 



A / 7"^ contrast between oil and water of 1.5 does not seem large. However, 
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this ratio appears as an exponent in an exponential function that describes the echo 
train amplitudes ( see e.g equation (5)). Consequently, ^^1^2 ^^^^^^^^ ^ 
provides an easily detectable amplitude difference between two fluids with identical T 
times. 

2 

[0052] The / ^2 ^^"^'"'^^^ exploited using several NMR methods, for example, 

saturation recovery (SR), inversion recovery (IR), or steady state method in 
combination with conventional CPMC pulse sequences. Among various methods for T 
^ evaluation, IR (which monitors the growth of the signal magnitude from 1 to +1) has 
more-dynamic range for signal detection than SR does (which monitors the growth of 
the signal magnitude from 0 to +1) (see FIG. 4A). In FIG. 4A. curve (a) depicts the 
exponential growth of an NMR signal magnitude after the saturation power has been 
turned off, and curve (b) depicts the exponential growth of an NMR signal magnitude 
after a 1 80 degree Inversion pulse. It is apparent from FIG. 4A that curve (b) would ' 
provide at least two fold better sensitivity In monitoring the exponential growth than 
would curve (a). As a practical matter, the measurements cannot start immediately 
after a saturation pulse is turned off. A delay time Is required for instrument circuitry 
to equilibrate before the detection can begin. This delay further degrades the 
sensitivity of the SR approach because most dramatic magnitude changes are lost 
during the delay time. For this reason, IR typically produces data that are about 3 fold 
better quality than those from SR. Therefore, the following discussion will use IR as an 
example. However, one skilled in the art would appreciate that the following 
discussion applies equally well to all other methods that can produce T ^ / 7"^ ratios, 
including SR. 

[0053] Using inversion recovery (IR) as an example, a pulse sequence for measuring the T 

I T ratio would be: (WT 1 80 RT CPMG) , , as illustrated in FIG. 5A, where WT is 
1 2 n 

the wait time (also called "polarization time") for spins to relax back to the steady 
state, RT is the recovery delay time during which the signals exponentially grow back 
to the steady state magnitudes, and CPMG is the typical CPMG pulse train (i.e., a series 
of 1 80-degree pulses with a fixed delay time, TE, between them). This pulse sequence 
will be referred to as an "IR-CPMG" sequence. 

[0054] 

In FIG. 5A, the first 1 80-degree pulse inverts the magnetization (i.e., onto the Z 
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axis). After RTde\ay the signal may be detected with a 90-degree pulse. As shown in 
FIG. 5A, the signals are actually detected with a CPMG sequence as a train of spin 
echoes. The detected spin echo magnitudes will be influenced not only by the 
relaxation process, but also the relaxation process, which is encoded by the 
various /?rdelay times. The IR-CPMG sequence as shown in FIG. 5A would detect 
signal components that have relaxed back to. the Z axis after RTdelay. Herein, the "IR- 
CPMG" would be used to refer to the sequence shown in FIG. 5A and any variant 
sequence which substantially include this sequence and achieve the same purpose. 
FIG. 5B shows a CPMG sequence coupled with a saturation recovery (SR) sequence, 
instead of an inversion recovery (IR) sequence. This sequence would be referred to as 
an "SR-CPMG" sequence. The SR~CPMG sequence as shown in FIG. 5B would detect 
signal components that have relaxed back to the Z axis after ^rdelay. As stated 
above, the SR-CPMG may be used instead of the IR-CPMG in the embodiments of the 
invention. 

[0055] Ideally, WT in any pulse sequence should be long enough for all spins to relax 
back to the steady state (i.e., fully polarized by the static magnetic field) before the 
next pulse cycle. However, time economy often dictates that a shorter WT is used, 
which may result in incomplete magnetization being detected. If insufficient WT is 
used, polarization of the spins by the static magnetic field may not be complete. The 
detectable signal intensity under this condition is given by: 

where / denotes the detectable signal intensity when the spin is fully polarized by 
o 

the static magnetic field (or fully relaxed back to the steady state). 



[0056] 



The measured signal amplitude ( /I . for the /-th echo from the CPMG train) using 
an IR-CPMG sequence would be: 



00 



where A is the full signal amplitude after full repolarization along the static magnetic 
field, RT\s the inversion recovery time in the IR-CPMG sequence (see FIG. 5A), TFis 
the delay time between the 1 80-degree pulses in the CPMG sequence, T ^ and are 
the longitudinal and transverse relaxation times, respectively, and 51s defined as: 

5' = l+XB*|l-«"^ J, 
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where /£is the inversion efficiency. The inversion efficiency, /£, is 1 .0 under ideal 
condition. However, /£is around 0.8 for a typical well logging tool. 

[0057] The above formulation assumes that signal decay due to diffusion is negligible, 
and the measured spin echoes decay with true . In addition to the decay 
during the CPMG pulse train, the measured spin echo magnitudes are also influenced 
by the ^rdelay times in the IR-CPMC sequence ( 7"^ relaxation). Alternatively, if RT is 
kept constant, the signal magnitudes would be Influenced by relaxation. 

[0058] Referring to curve (b) in FIG. 4A, it is apparent that there exists a critical RT, for 
which the recovering signal would vanish and become undetectable (i.e., >4 . = 0). This 
^roccurs at where curve (b) intercepts the time axis in FIG. 4A. Because different 
fluids have different T ^ values, each fluid may have a different ^7"time when the 
signal would become undetectable. Referring now to FIG. 4B, if a fixed RT delay time, 
e.g., 7"^ , is used in the IR-CPMG sequence, some signals would have recovered 
sufficiently to produce a positive signal magnitude (curve (c)), some might become 
undetectable (curve (d)), and some may produce a negative magnitude (curve (e)). 

[0059] According to equation (1 1), >4 -OwhenT", = RT / /n ( S). Thus, T ^ - >?r/ 

I Ic Ic 

0.59 for 5 = 0.8, or T = RT / 0.69 for 5 = 1 .0. Signal components with T < T 

Ic y M 1 Ic 

(curve (c) in FIG. 4B) will appear with A ^ > 0, and signal components with ^ ^ ^ ^ 

(curve (e) in FIG. 4B) will appear with A < 0, while signal components with T = T^ 

i 1 Ic 

(curve (d) in FIG. 4B) will be undetectable. 

[0060] FIG. 4B illustrates the node where >1 . = 0 in T^ measurement. This node can also 
be located using the first echoes in the CPMG measurements using the IR-CPMG 
sequence shown in FIG. 5A having different ^rdelay times. In well logging, a train of 
CPMG echoes are typically recorded for each ^rdelay. The node in the measured data 
where the intensity is zero will depend not only on T^ but also on T^ . Therefore, the 
position of this node (where /4 . = 0) in the T^ spectrum is determined by the T ^ / T 
^ ratio. By recording echo trains for several RT, it is possible to invert the data for two 
components with two ^^1^2 ^^^^^^ neighborhood of this node. The simplest 

measurement suite would consist of two measurements: one with a standard burst 
CPMG without the inversion recovery sequence and the other with an IR-CPMG with RT 
< WT, 
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[0061] For data inversion, if the / ratio for water is known, inversion routines, 

such as the Window Processing method as disclosed in U.S. Patent No. 5,291 ,137 
issued to Freedman, for normal CPMG data suites may be used. Only the data model 
for the experimental data has been changed to that described above. Data inversion 

would yield T distributions for oil and water directly. From these T distributions, 

2 2 

porosity fractions (water-filled porosity, oll filled porosity, and total NMR porosity), 
volume fractions (bound fluid volume and free fluid volume), constituent oil 
viscosities, viscosity or diffusion constants of the oil mixture or brine, and 
permeabilities can then be calculated. Methods for deriving formation properties from 
the distribution are known in the art. For example, U.S. Patent 6,229,308 Bl 
provides details on how this might be achieved. 

[0062] FIG. 6 illustrates key elements in the embodiments of the present invention. First, 
a suite of (i.e., a plurality of sets of) NMR measurements are acquired (shown at SI ) 
with at least one parameter selected from wait time or polarization time (WT), recovery 
time (RT), inter-echo delay time (TE), and number of echoes (NE) varied between 
different sets of measurements, wherein each set of measurements includes 
information on the longitudinal ( T ^ ) and transverse ( T ^ ) relaxation process. For 
example, the measurements may be acquired with the IR-CPMG or SR-CPMG sequence 
(as shown in FIGs. 5A and 5B), or any variation of these sequences so long as the 
collected data include information on both T ^ and 7"^ . The parameter is preferably 
so varied that the T ^ and /or process affects each NMR data set to a different 
extent. The sets of NMR data may be acquired with any NMR logging tool, wireline or 
logging-while-drilling. In addition, these measurements may be made with core 
samples removed from the formations in the laboratory. 

[0063] Once NMR data are collected, a formation model is generated (shown at S2) for 
inversion analysis. This model may include oils, brine, OBMF, and/or gas. In addition, 
within each fluid type (e.g., oils), the model may include multiple components to 
account for the constituents in the mixture. The model may include estimated signal 
magnitude, longitudinal relaxation time, transverse relaxation time, ratio of 
longitudinal to transverse relaxation time, diffusion constant, and viscosity for each 
constituent. This model is then optimized (shown at S3) to fit the collected data. The 
optimization process typically includes calculating a synthesized NMR data set based 
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on the model; comparing the synthesized NMR data set with the NMR measurements; 
and adjusting the model and repeating the calculating and comparing processes, if 
the difference between the synthesized NMR data set and the NMR measurements is 
greater than "a minimum." One skilled in the art would appreciate that "a minimum" 
for optimizing a formation model as used herein may be a true minimum of the 
difference between the calculated values and the experimental values. Alternatively, "a 
minimum" for the purpose of formation optimization may be a preset value, below 
which the difference between the model and the formations is considered acceptable. 
The optimization process may use any algorithm known in the art, for example 
PoweH"s nonlinear Constrained Minimization Algorithm cited In U.S. Patent No. 
6,229,308 Bl , or Monte Carlo simulation. 

[0064] The optimized formation model Is then output (shown at 84) with the computed 
parameters for each constituent fluid. These parameters may include or 7"^ 
relaxation time, diffusion constants, and constituent viscosity. From these constituent 
parameters, the desired formation properties can then be calculated (process SS). As 
stated above, the formation properties may include: porosity fractions (water-filled 
porosity, oil-dilled porosity, and total NMR porosity), volume fractions (bound fluid 
volume and free fluid volume), constituent oil viscosities, viscosity or diffusion 
constants of the oil mixture or brine, and permeabilities. 

[0065] The utility of embodiments of the present invention will now be illustrated with 
the following examples. In which computations are performed with a software 

TM 

package such as that sold under the trade name of MatLab by The MathWorks, 
Natick, MA. 

[0066] 

First, two identical T ^ distributions for oil and water (representing the most 
difficult case to resolve) were created with 7"^ values in the range from 2.5 to 2 5 ms. 
The oil filled porosity is 1 0 pu, the water-filled porosity is 20 pu. The ^^1^2 ^^^^^ 
was set to 1 .0 for the oil component and 1 .5 for the water component. These 

distributions of T were observed with 3 suites of 5 IR-CPMG or SR-CPMG and 1 

2 

CPMG sequences, having the following parameters: 
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Suite 1 Long WTIR Suite 






.IE 


OS 


0,8 


0.8 


08 


08 




WT 


10.000 


10.000 


10,000 


10.000 


10,000 


10,000 


RT 


1 


■JlQ 


10 


-Jvao 


100 


NA* 


TE 


05 


0.5 


0.5 


0.5 


0.5 


05 


NE 


500 


500 


500 


500 


500 


500 



Parameter is not applicable, since this represents the normal CPMG sequence. 







Suite 2- LongWT SR Suite 






IE 


00 


00 


00 


00 


0.0 


NA 


WT 


10,000 


10,000 


10,000 


10,000 


10,000 


10.000 


RT 


1 


VTo" 


10 




100 


NA 


TE 


05 


0.5 


05 


0.5 


0.5 


0.5 


NE 


500 


500 


500 


500 


500 


500 








Suite 3 Opttmal WT IR Smte 






IE 


08 


08 


08 


0 8 


08 


NA 


WT 


10 


-,^00 


100 . 


jCooo -xjlo.ooo 


10.000 


RT 


1 


JlQ 


10 


>^oo' 


100 


NA 


TE 


0.5 


05 


05 


05 


05 


05 


NE 


500 


500 


500 


500 


500 


500 



Note that RT = 0 is used as a flag to describe a standard CPMG preceded by WT. All 
times are in milliseconds (ms). Each sequence was repeated only once with 500 
echoes collected (NE = 500). 

[0067] 

One pu of Gaussian noise was added to the experimental data. Then, 2 5 noise 
realizations were run for each suite. The data were then inverted using the nn/sO 
TM 

function in MatLab . This function provides an inversion with positivity constraint. 
The data were either inverted directly or with window sums. The window sums 

contained the following echoes: 1 , 2-3, 4-6, 7-1 0, 1 1 -1 6, 1 7-25, 26-39, 40-60, 61 - 

91, 92-138, 139-208, 209-314, 31 5-500. The following tables show the results from 

the inversions (using "direct"or "window sums" (WS) approach): 
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Oil 


Water 








Suite 


Tail (diiect) (WST) 




1 


73±22 


7011.3 


2 


7.312.4 


6.9 11.2 


3 


65116 


7 2 1 1 .0 



Poroaty 







Oil 


Water 






= 0 10 





0 20 


Suite 


4) (direct) 




4> (direct) 


4) (WS) 


1 


0 U±002 


0 U ±G 01 


0 201002 


020+ 002 


2 


0 n±0 03 


0 11 ±0 03 


0 2010 03 


0 191003 


3 


0 1010 02 


0.1010.01 


0201002 


0 2010 02 



"WS"indicates inversions performed with window sums. 

[0068] These results show that the IR-CPMG suites (suites 1 and 3) provide better quality 
data than does the SR-CPMG suite (suite 2), as evidenced by the smaller standard 
deviations. As stated earlier, inversion recovery (IR) experiments monitor the growth 
of signals from 1 to +1 , while the saturation recovery (SR) experiments monitor the 
growth of signals from 0 to +1 (see FIG. 6A). Thus, it Is expected that the IR-CPMG 
suites would provide better quality data than would the SR-CPMG suite. 

[0069] A comparison between results from suites 1 and 3 reveals that suite 3, which has 
optimized WT values, produces better results. This is probably due to the spread of 
WT values that provide more accurate fitting of data to the model. This is so because 
the amplitude of each echo includes the term 

which can be better fitted with WT values spanning the 7"^ value (i.e., from WT < T ^ 
to WT > T ^ ), as compared with a single WT value in suite 1 . In addition, the 
parameters in suite 3 permit faster signal acquisition (than suites 1 and 2). The faster 
acquisition will allow more signal averaging/stacking within the same period of time, 
hence better signal-to-noise ratio (SNR). 
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[0070] The above simulation has assumed ^ ^ ^ I ^ ^ ^^^^^ ^ ^^'^ ^ ^1 ^ ^2 

ratio of 1 .5 for water. To investigate the effects of T ^ / T ^ ratio mismatch on the 

results, simulations were run with various T ^ / T ^ (water) ratios, ranging from 1.5 to 

2.4 in steps of 0.1 . The same T ^ distributions for oil and water were used. For each 

ratio, 1 00 noise realizations were averaged- Results for porosities and T are 

21m 

shown in FIGs. 7A and 7B, respectively. 

[0071] As shown in FIG. 7A, the porosity results are relatively sensitive to the ^^1^2 

ratio. A change of T ^ j T ^ from 1.5 to 2.4 resulted in an about 2-fold change in the 

porosities. In contrast, FIG. 7B shows that T . is less sensitive to T IT changes. 

21m 12 

A change of / T ^ from 1 .5 to 2.4 resulted in less than 40% changes in the T . 

1 2 21m 

values. 

[0072] Because the results depend very sensitively on the T ^ / T ^ ratio, this quantity 
needs either to be known a priori or it needs to be determined out of the data set as 
well. In principle, the data fit should be best for the correct ^^1^2 
error function (square sum of misfit between measurements and reconstruction) 
should be at a minimum for this value. In case of a linear model, this function can be 

written as: f = I ( mdata ( A/x spec)) ^ , where mdata are the measurements, M 
describes the linear model, spec 'xs the vector of model parameters, and >i/x spec 

describes the reconstructed NMR data based on the model. The summation goes over 

all measured data. 

[0073] FIG. 8 shows that there is, indeed, such a minimum that can be used to fix the T 
^ I T ^ ratio of water. The results in FIG. 8 are obtained with inversion of the window 
sums as described above. As can be seen from this inversion (which is averaged over 
1 00 noise realizations), the minimum is very flat and thus not well defined. However, 
such minimum may be better defined by averaging over large intervals with similar 
lithology (i.e., similar surface relaxivity). 



[0074] 



Instead of using a common T ^ / T ^ ratio for each fluid averaged over a large 
depth interval as illustrated above, it Is also possible to determine the T ^ / 7"^ ratio 
for each depth level in the inversion using a non-linear fit routine, e.g., the inversion 
method of Freedman (G.B. Patent No. 2,338,068 or U.S. Patent No. 6,229,30861) or 
the above cited non-linear optimization routine by Powell. However the weak 
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dependence of the IR-CPMG data suite on the ^^1^2 "^^^'^ makes the accuracy of 
the determined value questionable. Thus, the inversion using the T ^ / 7"^ ratio for 
each depth level may not be more accurate than that using a common T ^ / T ^ ratio. 

[0075] The above discussion presumes a homogeneous field. In an inhomogeneous field, 

the echo train decay constant may have contributions from both T ^ and 7"^ , and 

these contributions may change with the echo number. Thus, the effective T ^ / T ^ 

ratio for water might be echo number dependent for logging tools. In addition, it is 

assumed that the ^^1^2 "^^^'^ constant for each fluid. For 7"^ values in the BFV 

(bound fluid volume) region, the T ^ / T ^ ratio is probably constant. However, this 

may not be true for all fluids. For example, for water in large vugs, surface relaxation 

will be unimportant and thus the ratio will approach 1 .0. If this situation Is suspected, 

it is advisable to use other types of measurements to confirm the results. In addition, 

if the TIT ratio changes with T , the function should be known or be 
12 2 

determined. This function could be determined from cores or from analysis of a 
similar non-reservoir interval (without oil). 

[0076] Advantages of the present invention include its applicability to spins with short T 
^ or 7"^ , where diffusion-based approach fails. In addition, with short T ^ and T ^ , 
the IR-CPMG sequence may be repeated quickly td gather good SNR data in a 
relatively short period of time. Furthermore, the invention is not limited by the types 
of NMR tools used for acquisition; it is applicable to a wireline NMR tool, a formation 
tester having an NMR module, or an LWD NMR tool. In addition, methods of the 
invention may be applied to NMR measurements obtained in a laboratory using a fluid 
sample, which may be any mixture of fluids whether formation fluids or not. 

[0077] The above discussion uses the T ^ / T ^ contrast as a parameter in the MRF 
method to separate fluids. Other parameters may be used in combination.* For 
example, diffusion contrast may be used together with the T ^ / T ^ contrast. These 
information may be included in one data set and simultaneously inverted for optimal 

for the answer products. If the diffusion is not taken into account, the echo spacing 
preferably should be short so that diffusion is negligible and the decay constant is the 
diffusion-free T ^ value. 
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[0078] While the invention has been described with respect to a limited number of 

embodiments, those skilled in the art, having benefit of this disclosure, will appreciate 
that other embodiments can be devised without departing from the scope of the 
invention as disclosed herein. Accordingly, the scope of the invention should be 
limited only by the attached claims. 
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